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Heavy-light physics with NRQCD 

UKQCD collaboration presented by C. T. H. Davies a 

a Department of Physics and Astronomy, University of Glasgow, 
Glasgow G12 8QQ, UK 

First results are obtained for B mesons using a heavy propagator calculated using NRQCD, and a light Wilson 
propagator. Results from 13 quenched configurations of size 16 3 x 48 at p = 6.0 give a value for fs of less than 
200 MeV and a B* —B splitting of 32(8) MeV. Superior signal/noise behaviour is observed over static propagators 
on the same configurations. No extrapolation to the b mass for the heavy quark is required. 



1. INTRODUCTION 

B physics is an exciting area of topical inter- 
est. Results for B decay constants and mixing pa- 
rameters are urgently needed by the experimental 
community and can best be provided by calcula- 
tions in lattice QCD. There has been some contro- 
versy, however, in interpreting the lattice results. 
The static approximation gives large values for fg 
with a large spread from group to group jjj . The 
poor signal/noise properties of this approxima- 
tion presumably give rise to this variation. Cal- 
culations to date using Wilson propagators for 
the heavy quark have required an extrapolation 
in the heavy quark mass to the b. The use of 
NRQCD for the heavy quark gives a result with 
good signal/noise properties right at the b mass 
and the method is presented here. 

NRQCD is an effective field theory, appropri- 
ate to heavy quark physics, which has had great 
success to date in the simulation of heavy-heavy 
systems, both bb and cc ||,[|. 

In that case the NRQCD action is viewed as an 
expansion in powers of v 2 /c 2 where v is a typical 
velocity inside the heavy-heavy system. Simu- 
lations have been done recently by the NRQCD 
collaboration using leading and next-to-leading 
terms in this expansion [^|j3|. The coefficients 
of the sub-leading terms can be taken at their 
tree-level values when the input gauge fields are 
transformed as 
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quette Q. What this transformation does is to 
include in each term the radiative tadpole correc- 
tions to all orders in perturbation theory. 

The success of the spectrum calculated §,§ 
gives us confidence in the action and in the co- 
efficients of each term. That the associated per- 
turbation theory is well behaved and understood 
is shown by a study of the zero and finite mo- 
mentum energies of the T as a function of bare 
quark mass Q . The bare quark mass that should 
appear in the action can then be fixed by setting 
the kinetic mass of a known state, say the T, to 
its experimental value. This gives a bare quark 
mass for the b, Mt,a = 1.7 at /3 = 6.0 ||. 

To study heavy-light systems we can use the 
same NRQCD action for the heavy quark as has 
been used in the heavy-heavy simulations ||. 
There are then no parameters to be fixed at all 
and all masses are obtained as predictions. 

For heavy-light systems the importance of dif- 
ferent terms in the action changes. The momen- 
tum scale inside the mesons is set now by the dy- 
namics of the light quark and, to leading order, is 
independent of M, the mass of the heavy quark. 
The importance of terms in the NRQCD action 
then depends only on their power of 1/M. Terms 
at 1/M 2 which were included at next-to-leading 
order in a heavy-heavy simulation are effectively 
suppressed. We can neglect them when working 
to 1/M for a heavy-light simulation ||. 

The results described here use an NRQCD ac- 
tion for the heavy quark of the following form 



mq is taken as the fourth root of the average pla- 
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The heavy quark propagators are calculated 
with the usual evolution equation and with 
fields transformed as in equation 1. The bare 
heavy quark mass was taken as 2.0 to simulate a b 
quark on lattices at /3 = 6.0. Subsequent work on 
the heavy-heavy spectrum as decribed above Q 
shows this value for the mass to be 15% high. 
This initial study used 13 16 3 x 48 lattices fixed 
to Coulomb gauge. 

The heavy quark propagators were combined 
with light Wilson propagators at two Ki values, 
0.154 and 0.155. k c is determined at this (3 value 
to be 0.157 ||. The light quark propagators were 
calculated from a <5 function source and given a 
\J2ki rescaling. 

The heavy quark propagators were calculated 
both from a S function source and from an ex- 
ponentially smeared source (exp(— r/ro), with ro 
= 3.0). There was one source per configuration. 
By combining either smeared or local heavy prop- 
agators with local light propagators we are able 
to study local-local, smeared-local and smeared- 
smeared correlation functions for the B meson. 
We also calculated heavy quark propagators in 
the static approximation (that is, using equa- 
tion 2 with all 1/M terms switched off) from the 
same smeared sources. This enables us to directly 
compare the static and NRQCD signal/noise be- 
haviour. 

2. RESULTS 
2.1. f B 

To make a B meson we used the naive lowest 
order vertex operator, i.e. we combined the 2- 
component heavy antiquark propagator with the 
2 upper spins of the light quark propagator. 1/M 
effects in the vertex operator were ignored at this 
stage. The zero-momentum B meson correlation 
function then has an exponential fall-off given by 
its energy in lattice units and an amplitude re- 
lated to the decay constant, Jb- For the NRQCD 
heavy quarks the signal/noise is such that it is 
possible to fit the local-local correlation function, 
see Figure 1 . This has previously been shown not 
to be possible for static heavy quarks 0. Fit- 



ting to a single exponential from time slices 12-48 
and using an svd algorithm with cut-off 1.0e-03 
to invert the covariance matrix, we obtain at ki 
= 0.155: 

E = 0.48(1), Z L = 0.11(2), 
X 2 /dof = 1.7/4, Q= 0.8. (3) 

Zl is the square root of the amplitude. V%Zl is 
equal to JbV^Tb in the f n = 132 MeV conven- 
tion, modulo a lattice-to-continuum renormalisa- 
tion constant. This renormalisation constant has 
only been calculated for an NRQCD action with 
no a.B term and with no uq transformation of 
the [||. A fairly large 0(g 2 ) contribution was 
found in that case, but a major component of 
it was the wave function renormalisation of the 
heavy quark. This is much reduced when tad- 
pole effects are taken care of using eq. 1 [|]|]] . It 
seems likely that the renormalistion constant for 
the action used here is then much closer to 1. 

Extrapolating linearly to k c for the light quark, 
we obtain 0.16(3) /Za for /b\/Hb in lattice units. 
This agrees with other UKQCD results in this 
region |l(J , for reasonable values of Za ■ 

The predicted value for fs is then less than 200 
MeV, using an inverse lattice spacing at (3 = 6.0 
of 2.0 GeV. In the quenched approximation the 
value for the lattice spacing depends on the scale 
at which it is being determined; the a -1 value 
from bb splittings being considerably larger than 
that from light hadrons The scales appropri- 
ate to heavy-light mesons are presumably much 
more like those of light hadrons, so we use an a -1 
from those calculations, as above. 

The static calculation requires the calculation 
of both smeared-local and smeared-smeared cor- 
relation functions. The smeared-local correla- 
tion function is much noisier than the NRQCD 
smeared-local function (see Figure 2), but the 
smeared-smeared correlation function is only 
slightly noisier than its NRQCD counterpart. 

The poor noise properties of the local static 
correlation functions are because of the absence 
of a mass in the heavy-heavy channel that ap- 
pears in the noise fTl}| , Smearing has the effect 
of producing such a mass by effectively introduc- 
ing the heavy quark potential at non-zero R into 
the heavy-heavy cross-terms. This is a wholly 
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Figure 1. The effective mass for the local- local 
correlation function for the B meson using an 
NRQCD propagator for the heavy quark. 



desirable effect for the static smeared-local cor- 
relation function and the noise is much reduced. 
For NRQCD the non-locality of the source has the 
effect of increasing the noise over the local-local 
function, although an earlier plateau is seen. In 
both cases smearing at the sink introduces more 
noise than having a local sink. 

We can fit a single exponential to the ratio of 
the square of the smeared-local to the smeared- 
smeared correlation functions for the static case. 
The amplitude is then the effective local-local am- 
plitude. We obtain at ki = 0.155 

E = 0.47(1), Z L = 0.21(1), 
X 2 /dof = 6.7/13, Q = 0.9. (4) 

Of necessity there was a more restricted range 
in time (15 - 25) than for the NRQCD fit. We 
believe that this accounts for the fact that the 
errors are not significantly worse. 

The value for /%/M is 0.27(2)/ Z A after linear 
extrapolation to the chiral limit in ki. This is in 
agreement with other UKQCD static results at 
(3 = 6.0 [ fL0|| , obtained using the improved clover 
action for the light quarks, although the Za val- 
ues will be slightly different (by 10%) in the two 
cases. The value for fg is around 300 MeV, much 



larger than for NRQCD heavy qaurks. 

2.2. M B 

Because of the form of the NRQCD action, the 
energy measured for a zero-momentum correla- 
tion function must be adjusted to give the mass 
of the meson. This adjustment by terms which 
are perturbatively calculable works very well for 
heavy-heavy mesons . For the B meson we have 

M B = Z M M b ~E + E. (5) 

All quantities are in lattice units. Zm is the heavy 
quark mass renormalistion and Eq is the energy 
shift. Using E Q = 0.22 and Z M = 1-14 @ gives 
Mb = 2.46(5) in lattice units after extrapolation 
to the light quark chiral limit. Multiplication by 
a -1 is the major source of error, giving a result of 
4.9(4) GeV, certainly consistent with experiment 
(5.3GeV). 

Another determination of Mb comes from the 
finite momentum propagators. A kinetic mass 
can be extracted from the ratio of the ener- 
gies of lowest non-zero momentum to the zero- 
momentum. With only 13 configurations the 
results are too noisy to give a value from this 
method. 

2.3. B* - B 

From the ratio of the local-local NRQCD cor- 
relators for the B* and the B we can extract a 
value for the splitting between B* and B. A sin- 
gle exponential fit at k; = 0.155 to the ratio from 
time slices 5-20 gives 

E = 0.016(3), A = 1.02(2), 
X 2 /dof = 1.9/3, Q = 0.6. (6) 

Both values of ki give the same splitting so no chi- 
ral extraploation can be done. Taking the value 
above and a -1 = 2.0(1) GeV gives a splitting of 
32(8) MeV, rather lower than the experimental 
value of 46 MeV. This may be a reflection of the 
rather low quark mass that we have used. If the 
splitting is proportional to I /Mb, then a quark 
mass of 1.7 rather than 2 would give a splitting 
of 38(8) MeV, consistent with experiment. Val- 
ues obtained using heavy Wilson fermions have 
been much lower than this 112]. 
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Figure 2. A comparison of the effective masses 
for the smeared-local correlation function for a B 
meson using NRQCD propagators (crosses) and 
static propagators (circles) for the heavy quark. 



2.4. A fc 

We have calculated the NRQCD smeared(h)- 
local(l) correlation functions for the The ef- 
fective mass plots are rather noisy, and the mass 
depends quite strongly on m. Chiral extrapola- 
tion gives a splitting between the Af, and the B 
of 0.28(8) in lattice units. Converting to physical 
units gives the rather imprecise value at present 
of 60(20) MeV. The experimental result is 36(4) 
McV. 

3. CONCLUSIONS 

The use of NRQCD for heavy quark propaga- 
tors can give useful results for heavy-light mesons. 
Once the terms in the action are fixed from the 
very precise heavy-heavy calculations that are 
possible, there are no free parameters left. Corre- 
lation functions for the B meson can be calculated 
directly at the b mass (fixed from T spectroscopy) 
and they have good signal/noise properties. 

This preliminary calculation on a small number 
of configurations at (i — 6.0 finds a value for Jb 
consistent with others using heavy Wilson prop- 
agators. This confirms from another source the 



existent of relatively large 1/M terms in j\[M at 
the the B. The usefulness of the approach here is 
that a direct comparison with static results cal- 
culated using the same smearing functions on the 
same configurations is possible. 

A value for Mb can be calculated since all pa- 
rameters in the action are already fixed. This 
agrees with experiment. In addition the value for 
the B* — B splitting looks promising. It is reason- 
ably close to experiment unlike values obtained 
up to now using heavy Wilson fermions. 
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